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We treat the problem of normally ordering expressions involving the standard boson 
operators a, a) where [a, a" 1 "] — 1. We show that a simple product formula for formal power 
series — essentially an extension of the Taylor expansion — leads to a double exponential 
formula which enables a powerful graphical description of the generating functions of the 
combinatorial sequences associated with such functions — in essence, a combinatorial 
field theory. We apply these techniques to some examples related to specific physical 
Hamiltonians. 
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1 Introduction 

The normally ordered form of an expression involving ladder or, more generally, 
field operators is defined as one in which all annihilation operators are moved to the 
right using the appropriate commutation rules. Such expressions play an important 
role wherever the Fock space representation is used. In field theory and in many 
body quantum mechanics, normal order enters the formalism through the Wick the- 
orem [1]. This enables us to represent field-theoretical m-point Green functions or 
statistical m-point correlation functions (both defined as the vacuum mean values 
of products of the field operators) by expressions within which all initially taken 
operators are contracted. Another advantage of the normally ordered operators is 
seen if their matrix elements are calculated in the coherent states representation [2] . 
If the coherent states are introduced as eigenstates of the annihilation operators 
such matrix elements automatically become functions of the complex variables and 
provide us the Fock-Bargmann representation of quantum mechanical quantities. 
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Quantum optics is one branch of quantum physics where this approach is a basic 
tool widely used in numerous applications. 

The intersection between the boson normal ordering problem and combinatorics 
was discovered more than thirty years ago. The first seminal result [4] was that, for 
both bosons and fermions, ordering a general string (i.e., a product of nonnegative 
integer powers) of the canonical creation and annihilation operators results in an 
expansion in which the coefficients are combinatorial numbers called rook numbers. 
The second seminal result [5] expressed the boson string (aJa) n , where [a, at] = 1, 
as 

n 

(o)a) n = S jTS{n,k)a) k a k , (1) 

fe=i 

that is in terms of classical combinatorial numbers known as Stirling numbers of the 
second kind. Both results, especially the second one, being the "physicists version" 
of the relation 

= £S(n,fc)^ (2) 

known to mathematicians much earlier, have inspired investigations which lead one 
to conclude that the relation between the normal ordering problem and combina- 
torics is not accidental. Generalizations of (1) to the normally ordered expressions 
being more complicated than (a^a) , [6] -[9], have led to the introduction of new 
classes of combinatorial numbers, whose properties were elucidated using standard 
methods of combinatorial analysis — recurrences, generating functions and graph 
representations. Normal ordering for noncanonical operators has also been inves- 
tigated, in particular g-bosons, leading to ^-generalizations of Stirling and Bell 
numbers [10], [11], [12]. This work has also revealed the relation between matrix 
elements taken in the coherent states representation and the Bell numbers, [13], 
[14], [6], [7]. 

As noted above, all these results indicate a fundamental connection between 
the general normal ordering problem and combinatorics. Most of field theoretical 
calculations employ operator expressions reduced to the normally ordered form. 
This leads one to believe that the methods of combinatorial analysis will be useful 
in understanding and solving problems of quantum physics. Our main goal in this 
work is to justify this statement. We shall provide the reader with a general method 
for constructing normally ordered expressions and shall explain how to link them 
to well-known combinatorial problems. We shall also present analogies between 
our methods and those of standard field theory, in particular Feynman diagrams. 
Finally, we will illustrate this approach using examples arising in one mode boson 
normal ordering. 

2 The product formula 

oo oo 

Let f(x) = fnX n /n\ and g(x) = 9nX n /n\ be two formal power series, 

71=0 71=0 

also called the exponential generating functions (cgf) of sequences {/„}^L and 
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{ffn}~ =0 : respectively. Then 



x=0 



x Tx ) 



x=0 n=0 



which is a straightforward consequence of the elementary relation 



d" /x m 
dx n \ ml 



(3) 



(4) 



x=0 



reducing the usual Cauchy product of series to the point-wise Hadamard product 
If we apply this result to a function F(w) of an operator w then for any inde- 
terminate A 

F(At2>) = F ( X -$-) e x ™ 



da; 

On taking the normal form Af[ ■ } of the both sides 



x=0 



Af[F{Xw)} = F ( A— ) Af(c~ 



(5) 



(6) 



x=0 



We emphasize that on the left-hand side above the functional and operator aspects 
are mixed while on the right-hand side they are distinct. The functional aspects 
are given by a (formal) series in usual derivatives while the operator aspects are 
described by an universal expression — namely the normally ordered exponential 
of w, in general a word w(a,a^) in terms of the operators a, a). This means that 
(6), which we shall call the product formula, enables us to reformulate the general 
normal ordering problem into a normal ordering of e xw . To calculate it explicitly 
still remains a non-trivial mathematical task but the problem is tractable for a 
large class of physically interesting examples. In this note we shall consider the 
cases where w(a, a*) is either a product of positive powers of a and a\ or a power 
of a + a) . 



3 The double— exponential formula 

Using the notation 

Af(e^) = :G^(x,a,a^):, (7) 

where the symbol : : denotes that the function G^(x,a, a') is normally ordered 
assuming that and a commute, [2], [3], we can rewrite (6) as 



Af[F(\i 



F ( A— ) :Gv,{x,a,o}): 



x=0 



For many physical applications the general form of F is given by 



F( X ) = CXp I J2 
\m=l 



X" 



(8) 



(9) 
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and G may also be written in exponential form 

/ oo 

:G^{x,a,a)) := : exp ^ V^\a, a)) 



751 



(10) 



Substituting (9) and (10) in (8) we arrive at the so-called double- exponential for- 
mula 



Af[F{\i 



(11) 



x=0 



which is a convenient starting point for our analysis. Eqn. (11) means that using 
only simple derivations we can find solution to the normal ordering problem. Thus 
the function F defines the sequence {L m }^ =0 while the operator w defines the 
sequence {V n }%L . Examples of the solutions to the latter problem are [9] 



w = a^a . 



w = {a^) r a 



Af[exp(xa^a)] = : exp [a^a(e x - 1)] : , 
Vn W \a, at) = at a , 
Af [exp (a; (at ) r a)] = 



=: exp 



a) a 



]T( a t)(''-i)>_iy- 



x T{n+\/{r-l))x« 



w = a + a 



t 



r(l/(r-l)) n\ 
yW(a,at) = ata(at)(«-^(r - l) ^^/^^ > 

^•^(o+o+JJ — . e a; 2 /2 e a;(a+a t ) . 

F 1 W (a,at)=a + at, V/ 2 (ti) (a, at) = i , 
vJt*\a,at) =0 for n > 2. 



(12) 



4 Multivariate Bell polynomials 

Applying the double exponential formula requires some effort. It is easy to see 
that it works effectively if we deal with monomials in both exponentials but leads 
to rather tedious calculations in more complicated cases. It is more practical to 
expand both exponents as formal power series; a general method for this is the 
theory of multivariate Bell polynomials [15], [16]. 

Multivariate Bell polynomials arose from the question of constructing the Taylor- 

~ x n 

Maclaurin expansion of the composite function f(g(x)). For any f(x) = /« — 

n=l 

and g(x) — } g n — - given as formal power series one gets 
; nl 

n=l 

°° n 

f(9(x)) = [fog}(x) = J2F n [f;g}^, (13) 
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where 

n 

F n [f;g] = ^B nk ( gi ,92,---,g n -k+i)fk- (14) 
k=l 

The coefficients B nk are certain polynomials in the Taylor coefficients gi called 
multivariate Bell polynomials, or sometimes more simply but unprecisely 1 ), Bell 
polynomials. 

The multivariate Bell polynomials are closely related to combinatorial numbers. 
They satisfy 

B n , k (g u g n - k+1 ) = E"^^ 9? 9? ■ ■ ■ 9^1 , (15) 

m n w-wy>] 

3=1 

where the summation ^ " is over all possible non-negative {f,} which are partitions 
of an integer n into sum of k integers, e.L, 

n n 

^jvj=n, ^2uj = k. (16) 

i=i 3=1 

From (15) and (16) one may show that the multivariate Bell polynomials satisfy, 
for a and b arbitrary constants, the homogeneity relation 

B n ^ k (ab 1 gi,ab 2 g 2 , ab n - k+1 g n _ k+1 ) = a k b n B n<k {g u g 2 ..., g n -k+i) ■ (17) 

Recalling (12) one expects that the multivariate Bell polynomials of especial use to 
us are those from exponential generating functions 

exp(g(x)) = 22Y n [g]— , 

"=! "" (18) 

Y n[g] = ^2B nk (g 1 ,g 2 , . . .,g n - k +i) , 
fe=i 

obtained from (13) and (14) for f k = 1, k = 1, 2, 

It may also be seen from (15) and (16) that i?„,fc(l, . . . , 1) are the Stirling 
numbers of the second kind. In such a case (14) gives 

exp(u(exp *-!))=£ £ S{n, k)u k ) , (19) 

n=l Vfe=l / ^ 

n 

1 ) The Bell, or exponential, polynomials are defined as B„(u) = Yl S{n,k)u k , i.e., as polyno- 

fc=0 

mials which coefficients are the Stirling numbers of the second kind and B„(l) = B n are the Bell 
numbers. 
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and 



oo 



exp(expx — 1)) = ^ B n 



(20) 



with B n denoting the Bell numbers. 



71=1 

M 



x 



The second useful case is exp yg\x + qm -^jy ) , i.e., the case when only two g's 
do not vanish. We have 



/ x M \ °° x n 

exp [gix + g M j^J = ^2 H n M) (9i,9M)—f 



n=0 

n l M \ „n~Mr-r 

\ ^ 9\ 9m 



(21) 



Hi M \ 9l ,9M)-n\ ^ {n _ Ml . mm r 



where Hn M \gi, gui) are called the two variable Hermite-Kampe de Feriet polyno- 
mials. They are generalizations of the standard Hermite polynomials, [17]. 

Another important property is the inversion formula. This states that the 
following two expressions 



Y n [g] = ^-Bn,k(ffij . . .,g n _ k+ i) 



(22) 



fc=i 



and 



9n = £(-l) J '(j - • • -lYn-j+M) , (23) 

j=l 

arc inverse to each other if this notion is understood in the following way: an 
arbitrary series {^[g]}^! can be obtained in the form (22) if we choose the series 
{gn}%Li only as in (23). It means that we are really able to change the double 
exponential formula into power series and vice versa. Many other properties of the 
multivariate Bell polynomials are also known, together with their explicit forms for 
some basic (elementary) functions. 



5 The coherent states representation and the analogy to field theory 

As mentioned in the Introduction using the coherent states representation of the 
normally ordered strings of a and allows us to dispense with operators and deal 
with functions of a complex variable z and its conjugate z* . Taking the coherent 
states mean value of the double exponential formula we get 



(z\Af[F(\w)] \z) = exp ( £ % A- • exp £ ^ (,,*•) *- 



rn—1 



, (24) 



x=0 



(w) 



which, for the case z = 1 and for V, 



(l|^[F(A^)]|l)=ex P (^^A 



Vn W \l, 1), becomes 



dx r 



exp 5> ( *^ 



(25) 
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identical with the counting formula of [18], used there in order to enumerate the 
Feynman-type diagrams in zero-dimensional analogues of the field theoretical mod- 
els. 

The field theoretical analogy may be pushed further on recalling the functional 
formalism of the field theory, [19]. A basic quantity which defines any field the- 
ory model is the generating functional of the Green functions. Physically it is 
interpreted as the vacuum-vacuum transition amplitude of the time-ordered expo- 
nential of the quantum field operator in the Heisenberg picture 4>h(x) coupled to 
an external current J{x) 



G(J) = 



Texp^i J d D xj) H (x)J{x) 



(26) 



from which the m-point Green functions are got as the ra-th functional derivatives 
with respect to J. Passing to the interaction picture one obtains 

G(J) = (0|Texp^i J d D x (<Si n t(<Mx)) + 0i(z)J(z))) |0), (27) 

where <Si n t is an interaction Lagrangian density. Rewriting G(J) in terms of a 
functional integral 



G{J) 



= J [#] 



exp (i (S + S-u* + <f>J)) , 



(28) 



with So = J d D x£o and Si n t = / d D x£i nt denoting free (bilinear) and interaction 
action functional for the field 4>(x), and 4>J = J d D x(j)(x) J(x), respectively, one 
notes that G(J) is an analogue of the partition function of statistical mechanics. 
Another expressions for the generating functional of the Green functions are those 
given using functional differential operators 



G(J) :: OX P ( lA A A) 0X . p(i [«?.„ („)+„./]) 











exp 




• exp 


--JAJ 
2 



(29) 



where A = A(x, y) is the causal Green function of the free field equation generated 
by a free action, S- mt is, as previously, an interaction action and we use abbreviations 



d D xd D y T ^ 7 - K A{x,y)- d 



S(j)(x) 5(j)(y) 
JAJ = J d D xd D yJ(x)A(x,y)J(y) , 

Jcj) = J d D xJ(x)<j)(x) . 



(30) 
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Equivalence of both formulas of (29) comes from the identity 

D T{J) = T W)«pM] 



(31) 

4=0 



satisfied by arbitrary functional D and T having formal Taylor expansions around 
zero [20]. Here we remark that while in the derivation of the first equality of (29) 
the normal ordering is extensively used, the second equality, as well as (31), may be 
obtained by manipulating the functional integral representation of the Green func- 
tions generating functional. Although much simpler that the standard advanced 
field theoretical methods, our approach give essentially the same formulae, which 
are directly applicable to the time evolution operator or to the partition function 
integrand [21]. This we consider as a strong argument in favour of the present 
method. 



6 The counting formula, graphs and combinatorial field theory 

In the effort to understand the meaning of perturbation expansions in quantum 
physics, both in quantum mechanics and in quantum field theory, it is important 
to know their large order behaviour. Solving the problem for the coupling constant 
perturbation series one finds that the number of the Feynman diagrams, contribut- 
ing order by order to the perturbation scries coefficients, grows factorially, [22], 
[23], [24]. The factorial growth occurs because of the combinatorial reasons which 
(at least qualitatively and for the large orders asymptotics) explain why such a 
behaviour does not essentially depend on details of the model if one disregards the 
problem of renormalons 2 ) . Enumeration of diagrams for more complicated models 
is echoed by that of their zero-dimensional analogues, [25], [26]. Taking all dia- 
grams of a given order to have equal values shows that this is the number of graphs 
which determines the divergence of the (renormalized) perturbation series. Being 
able to estimate these numbers order by order we may investigate a perturbation 
series, estimate its asymptotic character and investigate resummation methods. 
The idea is to search for generating functions of given sequences of (combinatorial) 
numbers, i.e., to look for functions which enumerate the same weighted objects as 
the initial formal series do. If found in analytic form such functions may be treated 
as generalized sums of formal, i.e., divergent and so analytically meaningless, se- 
ries. We shall call such an approach a "combinatorial field theory" emphasizing 
that the notion should be understood in a two senses. The first sense is to obtain a 
generating function for the perturbation series, using any of the standard methods. 
The second, or inverse sense, is to construct interaction which leads to a postulated 
set of combinatorial numbers — following the statement of reference [18] that given 
a sequence of numbers {a n }, it is always possible to find a set of Feynman rules 
that reproduce that sequence. 

2 ) Renormalons arc a peculiar subclass of the Feynman diagrams which renormalized values 
grow factorially with the order of the graph. Renormalons appear in renormalizable models and 
are known to influence perturbation expansions of the nonabelian gauge theories. 
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The zero-dimensional analogue of (29) taken for J = (i.e., the generating 
functional for the bubble Feynman graphs in the field-theoretical jargon) is a par- 
ticular case of the right hand side of (25). So the latter formula may be considered 
as a recipe which counts the Feynman-like graphs in a model for which the poten- 
cy x n 

tial is 7 V n —r and for which we deal with " multileggcd propagators" of strenghts 
L m , [26], [27]. Using (19) and (20) we have for g's related both to {L m } and {V n } 



e u 9(x) = ^2 *J2u k Ynk([g]) = £ -Yndglu), (32) 

n=0 U - k=l n=0 U - 

and for (25) we get 

(l\Af[F(Xw))\l) = Z(L,V,X) = J2ML,V)^, (33) 

n=0 

using a notation 

A n (L,V) = Y n ([L])-Y n ([V}). (34) 

The standard example which illustrates how the counting formula works is the 
the zero-dimensional analogue of the anharmonic oscillator or (f> 4 model. It gives 

/ad 2 \ (-gx±\ l ^ T (2n+\) ( -gg\ n 



This is a consequence of 

- 1 T n = 



exp J xn = H "( x > y) ' ( 36 ) 

where H n (x, y) = Hn ' (x, y) arc the Hcrmitc-Kampe de Feriet polynomials, quoted 
in the Sect. 4, (21). For x = they take nonvanishing values only for even n 

= (37) 

The expansion in (35) may be equally considered either as the expansion in the 
number of vertices, i.e., in g n or as the expansion in the number of lines, i.e., in 
a". In (33) the expansion parameter is A which counts derivatives and, in a graph 
representation, produces lines or "propagators". We shall adopt this convention in 
what follows; however, each approach may be translated into the other using (31). 

The series at the right hand side of (35) is divergent which poses the question 
as to its meaning. The answer is given if one notices that the same series is got 
when one calculates a formal power series expansion in g for the integral 

'^^T^tiy^Ta^-i)- < 38) 
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The integral above does not define a function analytical for g = because it does 
not exist for g < 0. Being well defined for g > it admits an analytic continuation 
to g G {C 2 \ [0, — oo)} which is given by 

I(a, g) = JL (-) exp (-) K 1/A (-) , (39) 

where K is the Macdonald function, i.e., a Bessel function of the imaginary argu- 
ment. This allows to give a meaning to the series in (35) — it is an asymptotic 
expansion of (39) for g — > 0, g e {C 2 \ [0, — oo)}. So the series under consideration 
may be treated as summable in a generalized sense (nonuniquely, modulo functions 
obeying zeroth asymptotic expansions) to the properly defined function. 

As explained in the Sect. 4 we are able to find {A n }^L either constructing the 
multivariate Bell polynomials for sequences {L}^ =0 and {V'j^Lo or by identifying 
factors in the exponential formula with known exponential generating functions. 
Formula (33) also gives A n a combinatorial interpretation. Namely, the number 
A n is obtained as the number of all, connected and disconnected, graphs with the 
same number n of labelled lines. The Feynman-like graphs representation is how- 
ever different from the conventional one of field theory For the first thing, we 
classify graphs using the number of lines, not vertices. This does not matter for the 
simplest examples having monomials in both exponentials, but for more compli- 
cated normally ordered expressions our method leads to an alternative description. 
We refer the reader to [9] and [28] for details; here we list only the rules of the graph 
construction and define those multiplicity factors which are necessary in order to 
understand pictorial illustrations of examples in the next section. They are: 

— a line starts from a white dot, the origin, and ends at a black dot, the vertex, 

— we associate strengths Vk with each vertex receiving k lines and multipliers L m 
with a white dot which is the origin of m lines, 

— to count such graphs we calculate their multiplicity due to the labelling of lines 
and the factors L m and Vk- 



7 Examples 



We illustrate our approach using as examples the w — a^a and w — a) + a. We 
first consider the normal ordering problem 



cxp(Aat a +^ T ( a t a ) M 



(40) 



related to Hamiltonians \a) a + (a) a) M , in particular for M = 2 to the Kerr 
type Hamiltonian 



H = \a T n (l + ^aJa) . 



(41) 



In our previously introduced notation 

L\ = 1 , L M = 1 for some M > 1 , and L m = otherwise, (42) 
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and 



n 



( x M \ 

x)=cxp(x + —J; V^ = l for n=l,2,. 



(43) 



In order to get the Taylor-Maclaurin expansion of F(x) we recall (21) and expand 
(43) in terms of the two variable Hermite-Kampe de Feriet polynomials (x, y) 

^wrE^M)^ ( 44 ) 



From (12), (33) and (34) this yields 



(45) 



where B n are the Bell numbers. Combinatorially B n count all the partitions of an 
n-set and H^ M \l, j^) count partitions of an n-sct into singletons and M-tons. 



For M = 2 



V2J 



V2 



1, 2, 4, 10, 26, 76, 232, .... 



(46) 



are the involution numbers expressible using Hcrmitc polynomials H n {x) and the 
initial terms of A n arc: 1, 4, 20, 150, 1352, 15428, .... For M = 3 the coefficients 
A n are: 1, 2, 10, 75, 527, 6293, . . . , etc. In both cases the series diverge, as may 
seen from the n — > oo asymptotics, [29] 



B„ 



exp(expr(n) — l) 
"(n)]" +1 ^/27rcxpr(n) 



(47) 



where r(n) ~ logn — log(logn), applying d'Alembert criterion to A n /n\. 

We now obtain the generating function of the sequence{A„/n!}^L . First con- 
sider the partition function integrand corresponding to the Hamiltonian related to 
(40) for the particular value z = 1: 



cxp 



x -(a^a) M 



M! 



CO » n 

= E^ M) (-/^)-<r 

■P (x\ + 
-13 (kX + 



e 



(x\) n 



dx exp 
1 1 



A I 



e^kl 

k=0 



(xX)_ 
Ml 

(k\) M 
Ml 



X 

\k=0 
f oo 

E 

Vfc=0 



5{x - k) 
kl 

(5(x-fc) N 



(48) 
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n=1: 



n=2: 



n=3: 



1 = 1 



1 + 1 



< 



< 



= 4 



= 20 



Fig. 1. Lowest order Feynman-type graphs for the Kerr-type Hamiltonian with n = 1, 2, 3 
lines. The number below each graph is its multiplicity. 



where we have expressed the Bell numbers as the Stieltjes moments of the so- 
called Dirac comb, equivalent to the representation of the Bell numbers by the 
Dobihski formula. Applying formal manipulations we obtain a series which is con- 
vergent for positive (3 and which is identical with the expression which we would 

get if calculated ^1 |A/" exp ^— f3{\a)a + ^jr{ a ^ a ) M )j] from the definition of the 
Glauber-Klauder-Sudarshan coherent state \z = 1). 

Calculation can be repeated for an arbitrary coherent state using generalized 
Dobhiski formulae, [31], for the Bell polynomials 



k />oo 00 

B fe (|z| 2 ) = ]TS(fc,0M 2i = / dxx k e~^ 2 ^2 
i=o •'° 1=0 



It leads to 



exp ( -/3(Aa+a + — { (a^a) M 



OO 



\z\ 2l 6(x-l) 



^HM[-p,M Bn {\z\>)^ = 

n=0 ^ 



(49) 



A' 



00 U|2fc 



k=0 



exp 



(k\) 



M 



Ml 



(50) 

As mentioned the series (50) has a physical interpretation — for (3 = 1/kT it is 
the partition function integrand for the Hamiltonians related to (40). Moreover, 
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the series (50) may be integrated with respect to \z 2 \ term by term which leads 
to the partition function expressed in terms of the Jacobi theta functions and 
their generalizations [30] as expected. This we consider as a further argument 
which justifies the procedure undertaken and suggests that the method will also be 
effective in more complicated applications. 

The second example is the word w = a + and the normal ordering problem 
of Af (exp jjjy (Xw) M ) . Thus we consider 

F{x)=exp(j^), M =1,2,3, .... (51) 

Because 

TV (e x ™) = : Gn>{x, a, a f ) : = : e x2 /V (o+ot) : , (52) 

we have 

V^\a,a s )=a + (J , {a, a<) = 1 , V { ™ ] (a, a*) = for n > 2 . (53) 

In order to get the expansion of (52) we use the modified Hermite polynomials 

M*)=(-^)"^(^), exp (to (54) 
In general we get 



n = 



( \ M d M \ / x 2 

Z M (L, V, A) - exp ^— — J • exp j^2.x + — 

_ ^ ^Mn(2) ^ A MN 



z=0 



n! M! 

n=0 



(55) 



(Mn)! 

from which it is seen that the values of A n are Amu — 7 — p; r ^Mn(2) and zero 

(M!)™n! 

otherwise. The particular cases for which the generating functions are known in 
closed forms are: 



M=l, An = h n {2) = 1, 2, 5, 14, 43, 142, 499, 1850, ... , n = 0, 1, 2, . . . 

2 



Z 1 (L,y,A)=cxp(2A+^ 



(56) 
and 

M = 2, A 2n = 1,5, 129, 7485, 755265, 116338005, 

z i/ ^„(2) /A 2 \" 1 / 2A 2 \ (57) 



13 



A. Horzela et al. 



<; 

f « 4 

<; 
<; 



* 1 =129 



Fig. 2. Lowest order Feynman-type graphs for the superfluidity-type example with 
n = 2,4 lines. The number below each graph is (multiplicity) x FJ fc (vertex factor)Vfc = 

2<5fc,i + 5k,a). 



known as the Doetsch equality. Physically the latter example corresponds to a 
(special case of) single mode superfluidity-type Hamiltonian H ~ (a + a*) 2 ), [32], 
while mathematically it is the solution to the normal ordering of the exponential 
of the general su(l, f)-Lie algebra element, [33]. 

A closed expression has also been found recently, [34], for the case M = 3: 

n=0 v 7 

1 /\ 3 A 3 . A 6 \ _ (I 5 3A 6 \ 

= (1 - WW ^ T 2 ^ O U'6 ; ' ; 2(l-0A3)3 j' 

, (58) 

1 — \/i — 4A 3 

where 4>(X) = and 2-F0 is a formal series for (2,0)-generalized hy- 

pergeometric (MacRobert) function. In contrast to the M = 1 and M = 2 cases 
the series in (58) is divergent; however it is asymptotic and Borel summable to the 
function of well-defined analytic structure — the Airy function [35] . 



n=4: 




1 x 16 + 6*4 + 3x1 = 43 



Fig. 3. Feynman-type graphs for the case M (exp(^-A(a + a^) 3 )) with n — 4 lines. 
The case Af (exp(Jj(Au;) 4 )) is interesting because it corresponds to the c/i 4 /4! 
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interaction. We demonstrated in the Sect. 6 that for its zero-dimensional ana- 
logue the counting formula leads to a divergent series. Solving the problem ac- 
cording to our scheme we do not follow the rule (55) but exploit the knowledge of 

TV (exp given by Eqn. (57). Using (6), (36), (37) and (57) we have 



Af (exp |y (w) 4 ^ = exp 
and subsequently 



4! dx 2 



TV (exp ( 



cw 2 )) 



x=0 



(59) 



Af 



exp 



l>=exp(^ 



= E 



71=0 



2 

4! dx 2 



(1 - 2x 2 )V2 

y\ 4 



exp 



4x 2 



1 - 2x 2 



x=0 



2(4!) 



(60) 

The series is divergent for any y ^ and according to our best knowledge up to 
now no closed form of its generating function has been found. Investigating its 
relation to the example of the Sect. 6 we find 



Af 



exp 



= lim/z 

z^0\ 



Af 



exp 



= lim exp 

z^O 



lim exp 



v n=0 



______ 

4! dx 4 J 

______^| 

4! dx 4 J 

r(2n+l) 



z = 



exp 2xRe z + 



x=0 



E^^ 2) ((2Rez),l) 



x 
n\ 

4! 



n=0 
2„.\4 



a:=0 



(61) 

which generalized sum equals to (39) in view of (35). A similar technique may be 
used to find expansions of analogues of the Green functions and to investigate their 
summability by Borel or Pade methods but this important and interesting topic is 
out of the scope of this note. 



8 Conclusions 



The main thrust of this article has been to show the intimate relationship be- 
tween the normal ordering problem and combinatorics, both of series and of graphs. 
We have related the graphical approach to that of Feynman diagram in field theory. 
Basing our analysis on a simple product formula, we have applied these combina- 
torial methods to investigate some simple problems of field theory and many body 
quantum physics. The close connection between the normal-ordering problem of 
the operators appearing in the theory and the associated combinatorial algebra is 
evident. These combinatorial methods open the possibility of investigating solu- 
tions to a variety of field theoretical problems. Differing from, and more simple 
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than, the standard methods, this approach encourages us to believe that it may 
fruitfully be applied to more complex problems than those which we have treated 
for illustration in this article, and that it may well become a supplement to the 
standard methods of field theory, which are at present dominated by complex and 
functional analysis. 

We have benefited from the use of the EIS [36] in the course of this work. One of us (PB) 
wishes to thank the Polish Ministry of Scientific Research and Information Technology for 
support under Grant no: 1P03B 051 26. 
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